POINTWISE ESTIMATES FOR SPLINE GRAM MATRIX 
INVERSES OF ORDER THREE 



MARKUS PASSENBRUNNER 



Abstract. In this article we prove a pointwise estimate for inverses of Gram 
matrices corresponding to spline systems of order three. We use this inequality 
to obtain that the orthogonal projection Pa/ of / G L^[0, 1] onto the spline 
space of order 3 with mesh points A converges to / almost everywhere, provided 
the maximal mesh width of A tends to zero. 
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1. Introduction 

Let k E N and A = {{Ni^k, Nj^k))ij=i be the Gram matrix of the B-spline 
functions Aj^fc, 1 < * < M, of order k corresponding to an arbitrary partition 
A of the unit interval [0, 1]. For the exact definition of B-splines, see Section [2] 
Moreover, let be the orthogonal projection operator from L^[0,1] onto the 
space spanned by the spline functions {Ni^k)iii- 

Exponential decay inequalities (away from the main diagonal) for the matrix 
'■= translate into boundedness properties of P^. For instance, the 
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inequality 

(1.1) \kA< 

maxj I supp Ni^k \ , I supp Nj^k \ } 

for some constants C, < g < 1 (uniform in A), is equivalent to the uniform 
boundedness of as an operator from L°°[0, 1] to L°°[0, 1]. (cf. |CieOO] ) In fact, 
Sliadrin showed in his paper [ShaOT] that for each fixed k, this uniform bounded- 



ness of holds, implying (1.1). This additionally yields that for every partition 
A, the corresponding orthogonal spline system is a basis in C[0, 1] (and L^[0, 1] 
for 1 < p < oo). 

In order to show more properties concerning convergence of orthogonal spline 
series (like almost everywhere convergence or unconditional convergence), it is 



natural to improve ( 1.1 ) in a way that the desired convergence property is implied. 



For example, a uniform inequality of the form 
(1-2) \b.J < ^ 



maXjAj<i<ivi I suppiVi,fc| 

is used in |CK97] and |GK04j to prove convergence almost everywhere (for weak 
Lebesgue points) and unconditional convergence in for 1 < p < oo of general 
Franklin series (that is orthogonal piecewise linear splines) 



The aim of this paper is to prove the uniform inequality (1.2) for piecewise 
quadratic splines (A; = 3). As a first application of this result, we show the almost 
everywhere convergence of Pa/ to / for f E L^, provided the maximal mesh 



width of A tends to zero. The main ingredients of the proof of (1.2) are 

(1) The Sherman-Morrison formula applied to rank- 2- updates of matrices and 

(2) the total positivity of the matrix A. 

The organization of this article is as follows: In Section 2 we recall some 
preliminaries that include the Sherman-Morrison formula for inverse matrices, 
the definition of B-spline functions and some results concerning total positivity of 
B-spline Gram matrices. Section 3 contains the formulation of our main theorem. 



that is the exact form in which (1.2) holds. In Section 4, we use the Sherman- 



Morrison formula to derive an expression for iterated calculation of inverse Gram 



matrices of B-splines, which is the main starting point of the proof of (1.2). In 
Section 5, we apply our method of proof to obtain the known inequality (1.2) 
for piecewise linear splines. The purpose of this presentation is to illustrate the 
basic steps that are essentially the same as in more general cases. In Section 6, 



we derive our main result (1.2) for piecewise quadratic splines. Section 7 gives a 
few applications of this result, including almost everywhere convergence of Pa/ 
to / for / G L'^[0, 1] as the maximal mesh width converges to zero. 

2. Preliminaries 

2.1. The Sherman-Morrison formula. We have the following formula for the 
inverse of a rank m perturbation of a given matrix A. 

Theorem 2.1. |SM50l IWooSO] Let A be an invertihle n x n matrix and U,V 
n X m matrices. If the m x m matrix 

1 + V^A-^U 
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is invertible, then we have 

{A + UV^)-^ = - A-^U{1 + V^A''U)-^V^A-\ 

Remark 2.2. In applications of this formula, m is typically much smaller than 
n. We will apply it for the choice m = 2. 

2.2. B-splines. Suppose we are given a partition 

(2.1) A = (si)f^i, = si < S2 < •■■ < < Si = 1 

of the unit interval [0, 1] and an integer k > 1. We denote by Pfc(A) the space of 
piecewise polynomials of order k with respect to the partition A, i.e. 

(2 2) P'fc(^) •= {/ ^ • f\{s^,sr+i) is a polynomial of degree k - 1 

for each l<z<L — 1}. 

Furthermore, we define the spline space 

(2.3) Sfc(A,a):={/GP,.(A) : / G C"^^ (s,) for all 2 < j < L - 1}, 

where {<Jj)^Z2 is a sequence of integers (smoothness parameters) such that < 
(Tj < k — 2 for all j. Corresponding to the partition A, we define an extended 
sequence of points A where we allow each knot Sj to occur more than once. For 
this, let fi = (/ii, . . . , fii) with fj,i = fj,L = k and 1< fii <k — 1 for 2<i<L — 1 
be a multiplicity vector. We denote the £^-norm of a vector by | ■ |. Then we define 
A = A(A,yu) to be the unique sequence of points (ti)i=i that satisfies 

|{1 < < l/^l : = Sj}\ = fij for all 1 < j < L and 
ti < ^i+i for all 1 < i < — 1. 

So, A is the ordered sequence of points in A where each point Sj occurs /i^ times. 

Corresponding to the sequence A = (ti)l=i, we define the B-splines Ni^k of order 
k as 

(2.4) N,4x) = i-l)HU+k - U)[ti, . . . ,t,+fc](x - l<t<\fi\-k=:M, 

where [ti, . . . is the standard (forward) divided difference operator with re- 
spect to the points t,, . . . ,ti+k and is the truncated power function defined 

as 

M+ := X[o,oo)(m), := u''ul hi k > 1. 
We remark that A^j i is the characteristic function X[u,ti+i) of the interval [tj,ti+i). 
The relation between Sjt(A, cr) and the defined B-splines corresponding to the 
sequence of points A (A, /i) is given by the following well known theorem, see for 
instance |Sch81] . Theorem 4.9. 



Theorem 2.3. Let k eN, A as in (2.1) and a a sequence of smoothness param- 
eters as above. Define the sequence fi = {k,k — 1 — (r2, . . . ,k — 1 — ctl-i, k) and 
let Ni^k be the B-splines corresponding to the extended sequence of points A(A,/i) 
defined above. Then 

s^im{Ni^kYiLi = Sfc(A,a), 
i.e. the B-splines are a basis in the space of piecewise polynomials with certain 
prescribed smoothness at the mesh points. 
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We collect classical properties of the B-splines Ni ^ in the following 

Proposition 2.4. The B-splines Ni^k, k > 1, admit the following properties. 

(1) Ni^k{x) > for allx G R, 

(2) suppiVi^fe = [ti,ti+k], 

(3) Y:f=lN^,k{^) = '^for all xe [0,1], 

(4) the recursion formula for B-splines: 

Ni^k+l = T^i,k + 7 — 7 ^j+l,A:, 

ti+k ~ ti ti+k+l ~ ti+1 

(5) the L^-norm of Ni^k is given by 

\\Ni,k\\i = Ni^k{x) dx = h±!L_Jl_ 

In order to keep future expressions involving distances of points in A simple, 
we introduce the following notation. 



Definition 2.5. Let A; G N, A given as in (2.1), /i a multiplicity vector and 
A = A (A, /i) = be the extended sequence of points as above. Furthermore, 

let = for i < and tj = 1 for i > \fi\ + 1. Then we define 

(2.5) (m, n)j := {mn)j := tj+rn - tj+n 

for integer parameters m, n,j. We define 77 : {1, ... , M} x {1, . . . , M} — )■ ]R_|_ to 
be an arbitrary function that satisfies the following axioms: 

(1) ri{i,i) = {k,0)i for all 1 < i < M, 

(2) = r]{j,t) for all l<t,j<M, 

(3) vU, n) < vU, n + l)< r]{j, n) ^tZV+i for all 1 < n < M - 1, 1 < j < n. 
Example 2.6. (1) The function 

(2.6) Viihj) ■= max (fc,0)/ 

iAj<l<i\/j 

satisfies the axioms (l)-(3) of the above definition. Indeed, axioms (1) 
and (2) and the left hand inequality of axiom (3) are trivially true for 771. 
The right hand inequality of axiom (3) is an immediate consequence of 
the two obvious inequalities 

{k - 1, 0)„,+i < {k, 0)„+i, {k - 1, 0)„,+i < (A;, 0)„ < r/i(j, n) 

and the definition of rji. 
(2) The function 

(2.7) V2ihj)-=tiVj+k-tiAj 



satisfies axioms (1) - (3) of Definition 2.5 Axioms (1), (2) and the left hand 
inequality of axiom (3) are again immediate consequences of the definition. 
The right hand inequality of property (3) follows from the elementary 
inequality 

a ^ 1 

(a + 6) (a + c) ~ a + 6 + c 
for positive real numbers a, b, c. 
Clearly, r]i{i,j) < mihj) < li-jlviihj) for all l<i,j <M. 
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Remark 2.7. Since the values for m, n in the expressions (m, n)j are not greater 
than 6 in our cases, we typically write for instance (10)j to denote (1,0),. 

2.3. Total positivity. We denote by Qm,n the set of strictly increasing sequences 
of m integers from the set {1, . . . , n}. Let A be an n x ra-matrix. For a, /3 G Qm,n, 
we denote by A [a; /3] the submatrix of A consisting of the rows indexed by a and 
the columns indexed by /3. Furthermore we let a' (the complement of a) be the 
uniquely determined element of Qn-m,n that consists of all integers in {1, . . . , n} 
not occurring in a. In addition, we use the notation A{a] (3) := y4[a';/3']. This 
notation is taken from |And87] . 

Definition 2.8. Let A be an n x n-matrix. A is called totally positive, if 

(2.8) det A[a; /3] > 0, for a, /3 G 1 < m < n. 

The cofactor formula bij = {—ly^^ det A{j]i)/ det A for the inverse B = 
{bi,j)ij=i of the matrix A leads to 

Proposition 2.9. Inverses B = {bi,j) of totally positive matrices A = {aij) have 
the checkerboard property. This means that 

{-iy^'bij>0 for all t J. 

We introduce the notation {f,g) := Jq fg for the standard inner product of 
two functions on the unit interval [0, 1] with respect to Lebesgue measure. 

Remark 2.10. The Gram matrix of B-splines A = ((A^i,fc, -^j,A;))fj=i of arbitrary 
order k and arbitrary partition A is totally positive (see |dB68] ). This is a conse- 
quence of the so called basic composition formula (Equation (2.5), Chapter 1 in 
|Kar6 8]) and the fact that the kernel Ni^k{x), depending on the variables i and 
X, is totally positive (Theorem 4.1, Chapter 10 in |Kar68j ). Thus the inverse of 
A possesses the checkerboard property by the above proposition. 



3. The main result 

Our main result is about the inverse Gramian of piecewise quadratic splines. 

Main Theorem 3.1. Let k = 3, A be an arbitrary partition of the unit interval 
[0, 1] and fi a n ar bitrary multiplicity vector. Moreover, let rj be any function as 
in Definition 2.5. Then, the inverse B = (&jj)f^=]^ of the Gram matrix A = 
{{Ni^3, Nj^3))f'j^^ satisfies the bound 

(3.1) \bi,j\<Ci^^ forl<z,j<M, 

where < q < 1 and Ci > are constants independent ofi,j and A. They can 
be chosen as q = (87/100)^/2 Ci = C{l + ^j^), where C = 12g-2(|)^ 

Remark 3.2. It was proved by Shadrin |Sha01j that the L°°-norm of the L^- 
spline projector is bounded independently of the knot sequence for each /c G N. 
For this assertion there are a number of equivalent properties in |CieOOj . One of 
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them is that the inverse B = {bi^jij 
satisfies an inequahty of the form 



^^^^ of the Gram matrix A 



(3.2) 



for 1 < 2, j < M, 



supp A^i,fc| V I suppiVj-fcl 

where the constants c^, < < 1 depend only on the sphne order k. The result 
of Theorem 3.1 is thus an improvement of this inequality for A; = 3, since in 
Viihj) = niaxj/\j<Kivj I supp A^'i^fcl, the maximum ranges over all indices between 
i and j in contrast to the maximum in (3.2). For an analogous result in the case 
k = 2, cf. Remark [Kil and [KS89J. 



The first tool for the proof of Theorem |3.1| is an iteration formula for inverse 
matrices presented in Section |4j To illustrate the method of proof, we first show 
an analogous estimate for piecewise linear splines in Section |5j The next step is to 
consider piecewise quadratic splines, that solely have double points, which is done 
in Section 6.2 This limiting case also provides some insight in the general case of 
piecewise quadratic splines. Finally, Theorem |3.1 is a corollary to Theorem 6.17 



proved in Section 6.3 



4. An iteration formula for inverses of matrices 



Let A 



\,j)f\=i for some M G N be an invertible M x M- matrix and define 



Ar. 



[(^i,j>i,j=l 



for Kn < M 



to be the nxn matrix consisting of the topmost n rows and the leftmost n columns 



A 



n+l 



'nxn 
,n\T 



of A. Furthermore set i?„ := if An is invertible and let Bn = (6"j)",_i. Then 
An+i may be written as the sum of two matrices as follows. 

An Onxl V (^0, 

is the column vector (a„+i^i, . . . , a„+i,„)"^ and is the column 
vector {ai^n+i, ■ ■ ■ ,o,n,n+i)'^ ■ Tn is a rank 2 matrix that can be written as the 
product Tn = UnV^ where f/„ and Vn are (n + 1) x 2 matrices and defined as 



(4.1) 



where f " G 







Un 







Onxl 
1 



Vn 



Onx 1 
1 







Theorem 2.1 and the above decomposition (4.1) of An+i in two summands 
yield the following 

Corollary 4.1. Let 1 < n < M. Additionally, suppose that An,Bn and f",^" 



are defined as above and set v 



■u". If An is invertible, an+i,n+i 7^ and 



an+i,n+i — V BnU ^ we have the following formula for Bn+i = A 



n+l 



B. 
Oi 



n 

xn 



Onxl 





+ 



v'^BnU 



n+l ■ 

BnUv'^Bn -BnU 
-V^Bn 1 



Q-n+l.n+l 

Proof. We compute the matrix products occurring in the Sherman-Morrison for- 
mula step by step and conclude 



Sn^Un 



Bn 
Olxn 



Onx 1 

*n+l,n+l 



Onxl 
1 



BnU 




Onxl 

^n+l,n+l , 
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T/^Q-irr _ /^Oixn l\(BnU 0„xl \_( a„+i,n+i 



Tc-lrr \-l _ / 1 '^n+l,n+l 



(1 + v^s-'u^r' 



1 - a;|i,„+ii;^S„M V-^'^^n^ 1 



Define x := (1 — a„:{:;^ ^^^v^Bnu) ^. We get further 

°"r)(-A. """1''"") 



™ / ^ '^n+l,n+l'" 

^ ' -v'^BnV 1 



1 T 



\yixn "n+l,rt+l/ 



f^5n -a~],„,.v'^B„u 



' —B„UV^ Br, BnU 
^71+1, n+1 

If we insert these expressions in the Sherman-Morrison-formula, Theorem 2.1 
conclude 



Bn 0„xl \ , 1 (BnUV^Bn -BnU 



we 



Oixn a^+^ n+ij a^+i^n+i - v'^BnU \ -V Bn BnU 

Bn 0„xA ^ 1 [BnUv'^Bn -BnU 



This proves the corollary. □ 



We now employ this corollary in the cases where the matrix A is a symmetric 
tridiagonal or a symmetric 2-band matrix. 
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4.1. Tridiagonal matrices. Let A he a symmetric tridiagonal M x M-matrix. 



Br, Onxl \ j_ l,n „2 



Using Corollary 4.1| on A, we get for \ <n < M — \ 
where C„+i is given by 



n+l, 



-6^„a„,„+i if 1 < i < n, j = n + 1 
-&",n«n,n+l if 1 < j < n, i = n + 1 
1 if z = j = n + 1. 



4.2. 2-band matrices. Let A 

means that a,' 



a. 



be a symmetric 2-band matrix. This 



— for |z — j| > 2. Using Corollary 4.1 on A, we get for 

2 < n < M - 1 



(4.3) 

where 

(4.4) 



n+l 



nx 1 



Bn 
Olxn 



"n+l, n+l 



and C„+i is given by 



(4.5) (a+l)^,, 



^(0'n+l,n+l — ^n,n'^n,n+l ~ '^Ki-\,rfl'n-\,n+\0'n,n+\ 
"n-l,n-l"n-l,n+iy 

'-6"„a„,„+i - 6"„„ia„_i,„+i if 1 < i < n,j = n + 1, 
-6"„a„,„+i - 6^„_ia„_i,„+i if 1 < j < n, i = n + 1, 



(C*n+l)j,n+l(C'n+l)j,n+l 
1 



if 1 < z, j < n, 
if i = j = + 1. 



The following two lemmas are independent of the special form of matrices con- 
sidered in the next sections. We will use them in Section 16.31 to estimate inverses 
of Gram matrices corresponding to splines of order 3. 

Lemma 4.2. Let A he a symmetric ^-handed M x M-matrix such that = A~^ 
is checkerboard for 1 < n < M . Then, the inequality 
(4.6) 



bnXl,n+l — ^'^n+l.n+l bn,n^ri,n+l (o-n.n+l 



2ctn,n'J^n— l,n+l 
Q-n— l,n 



) 



, '2n,n+lQ'n— l,ra+l 
Q-n— l,n 



X -1 

"n-l.n+l^n-l.n-ll-"- ' "n,n"n-l,n-l'^n~l,n) J 



holds for2<n< M -1. 



Proof. By (4.4) and the checkerboard property of Bn, b'^+in+i is given by 
(4.7) 



Ln+l 
"^n+l,n+l 



(Cfn+l,n+l ^n,n'^n,n+l + „ |a„^„+ia„_i_„+i ^n-l,n-l'^n-l,n+l ) 



Another consequence of the identities (4.3) - ( [4^ ) is 



b"^' - Kj ^ + ^n,n(^"n-2'^n-2,n + bl^^^an-l,n) ib'^^n-2^n-2,n + &j,„iian-l,n) 
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for 1 < i, j < n — 1. The choice i=j=n — 1 in these equations together with 
the checkerboard property of Bn yield the estimate 

(4-8) ^n-l,n-l — ^n-l,n-l(l + ^n,n^n-l,n-l^n-l,n) ■ 

The defining property of the inverse matrix Bn = A^^, the fact that is 2- 
banded and the checkerboard property of i?„ again imply 



(4.9) 



l^ji-l,nl ~ ^n~l,n{^n,n^n,n + ^n-2,n^n-2,n 



1) 



Finally, inserting (4.8) and (4.9) in (4.7) yields the conclusion of the lemma. □ 



Remark 4.3. The crucial fact about the foregoing lemma is that inequality (4.6) 



depends only on the matrix A and on the rightmost and bottommost entries 
KiwKi~-in-i of the matrices Bn, Bn-i respectively. 



Lemma 4.4. Let A be as in Lemma \4^ and 2 < n < M . Then, if I < j < n — 1, 
the estimate 

(4.10) \b]n\ < \bTn'^,Knan-l,n 

holds. Additionally for 3 < n < M and l<j<n — 2we have 



(4.11) 



\bln\<\bTn-lKA^n-l,n- 



(^n—2,n(^n—l,n—l 



a. 



n—2,n—l 



Proof. First, we note that ( 4.3[ ) - (4.5) yield 



(4.12) bln = -blnibrn-2^n-2,n + b-' 



if 1 < J < n - 1. 



Since Bn is checkerboard, (4.10) follows for 1 < j < n — 1. By the defining 
property of B„ 



An we obtain 



(4.13) 



LTl-l 



S^^n— 3,n— 1 



+ b 



n-l 
j,n- 



2^n—2,n- 



-1 + bln\an~ 



l,n-l 



Thus we get, under the assumption 1 < j < n — 2 
(4.14) ' 



a„l2,n-l(«n-3,n-l|&",„i3l + a„-l,„,-l 1 | ) 



— ^n-2,n-l^n-l,n-l\bj^n-l 



since Bn-i is checkerboard. Now use the checkerboard property of Bn in (4.12) 



and insert estimate (4.14) in (4.12) to conclude the assertion of the lemma. □ 



5. PlECEWISE LINEAR SPLINES 



We now apply Corollary |4.1| to the case of piecewise linear continuous splines 
to get geometric estimates for the entries of their Gram matrix inverse. We let 
^ = {si)^=i be a partition of the unit interval [0, 1] as in (2.1) and set the k = 2. 



The multiplicity vector fi in the case for piecewise linear continuous splines is 
given by 

/x=(2,l,l,... ,1,1,2). 

We construct the extended sequence of points A = A(A,/i) = (tj)|=i as in Sec- 
tion 2.2 and note that in the underlying case, M := — k = L. Furthermore, 



we construct the B-spline functions A^j = A^j 2 for 1 < z < M as in Section 2.2 
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that is in this section Ni is the unique piecewise hnear continuous function on 
the unit interval [0, 1] such that 



6ij for 1 < j < -L. 



Now we consider the Gram matrix A = {cii,j)ij=i = {{^i, Nj))ij=i obtained 
from these functions. Using the special form of the piecewise linear B-splines Ni, 
we get 

(5.1) tti^i = (20)i/3 if 1 < 2 < M, 



(5.2) 
(5.3) 



a 



(21)i/6 




if 1 < i < M - 1, 
if |z - j| > 1. 



Note that A is symmetric and tridiagonal, therefore we may apply the formulas 
of Section 4.1 to A. If we do and insert the expressions (5.1) -(5.3) for the Gram 
matrix A, we obtain 

(5.4) = 

where 



(5.5) 



(a 



n+l)i,j 



fW,n(21)n/36 ifl<^,j<n, 

-6[;„(21)„/6 if l<z<n,j=n + l, 

-b]j21)n/6 iil<j<n,t = n + l, 

1 a i = j = n + 1. 



Lemma 5.1. Let the matrix A be as above and 1 < n < M. Then we have the 
estimates 

3 , 3 4 



(5.6) 



(20), 



< 



[(10)„ + (21)„ - (20), 



for the rightmost bottommost element b^^ of the matrix = A^^ . 

Proof. We proceed by induction on n. For n = 1, we have by definition of Bi 

3 



and (5.1 ) 



(20), 



6i,i=3(-(10)i + (21)i 



since (10)i = and thus equality on both sides of ( 5.6[ ). We know from Remark 
2.10| that A is totally positi ve, so a fortiori An is totally positive for every 1 <n < 



M . Therefore Proposition 



2.9 



yields > 0. So we see the lower estimate 

in (5.6) immediately by glancing at formula (5.4) for For the upper 

estimate, we obtain as a consequence of the inductive hypothesis 



tn+l 
"n+l,n+l 



= 3(^(20)„+i 
< 3f(20)„+i 



-3 



^(^(10)„ + (21)„)"\21)^ 
7(21)n 



-1 



(10)„+i + (21)„+i 
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□ 



Lemma 5.2. Let the matrix A be as above and I < n < M . Then, the elements 
of the rightmost column of satisfy the geometric estimate 



(5.7) 
where q = 2/3. 



16" I < 



for all 1 < j < n, 



Proof. We infer from (5.4) and (5.5) that 



"j,n+l 



Q^n+l,n+l"j 



n+1- 



Invoking Lemma 5.1 for b^'^\ we get further 
(5.8) ,/2(10)„+i 



un+l 



< 



,6" I 
3(20)„+i' 



Now we note that by Lemma 5.1, inequality (5.7) is true for n = j. For general 



n > j, inequality (5.7) follows from (5.8) by induction using axiom (3) of rj in 

□ 



Definition 12.51 

Theorem 5.3. Let the matrix A be as above and 1 < n < M. Then, the elements 
of Bn = A^^ satisfy the geometric estimate 

36 gl'-^l 



6". I < 

5 r]{i,j] 



for all 1 < i, j < n, 



where q = 2/3. 



Proof. We first observe that it is sufficient to prove the theorem for the parameter 
choices i < j < n — 1, since Bn is symmetric and the case j = n is already covered 



by Lemma 5.2 Equations (5.4) and (5.5) yield 



6". 



tn— 1 _|_ rn ^n— 1 1 
f^-- "T "n,n'^i,n-l'^j,n-l gg ■ 



SO by induction 



Using now Lemmas 5.1 and 5.2 we get 



36 



n-l 



< 



4g- 



4 



q'^^-M^,j)mii 

9t^^(^,0^(j,0(20)m 



n-l 



Since (10)i+i < (20)i+i, (10),+i < (20); < r]{j,l) and r]{i,j) < r]{i,l) for j <l< 
77, — 1, we estimate this from above by 
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proving the theorem. 



□ 



Remark 5.4. An exponential estimate related to the one of Theorem |5.3| for 
the inverse Gramian B = {bij)f^j^i of piecewise linear spline functions was first 
proved by |Cie66] . Another proof, using the B-splines Ni can be found in |KS89] 
yielding the estimate 

2-l*-il 



< 4, 



supp A^j I V I supp Nj I 



The same argument in fact yields 



Vi[hJ) 

which is better than the assertion of Theorem |5.3[ Nevertheless, this proof strongly 
uses the tridiagonal structure of the underlying matrix A and thus cannot be gen- 
eralized to Gram matrices of higher order splines. 

The reason why we presented Theorem |5.3 anyway, is that the basic structure 
of its proof is the same as in the more general cases of piecewise quadratic splines 
in Section [6l 

6. Piecewise quadratic splines 



In this section, we apply Corollary |4.1| to the case of piecewise quadratic splines 
to get geometric estimates for the entries of their Gram matrix inverse. We let 
^ = {si)f=i be a partition of the unit interval [0, 1] as in (2.1) and set k = 3. The 
multiplicity vector n in the case for piecewise quadratic continuously differentiable 
splines is given by 

/i=(3,l,l,... ,1,1,3). 

We construct the extended sequence of points A = A(A,yu) = and the 

B-spline functions Aj = A'j s for 1 < i < M = |/i| — 3 as in Section 2.2, that is 





U = U 



U <tA< 



t 



1. 



6.1. The Gram matrix. We now calculate the Gram matrix of these B-splines. 
There is a standard formula for the inner product of B-splines of arbitrary order 
k involving divided differences (see for instance jSchSl] , Theorem 4.25), but in 
this section we prefer to calculate the Gram matrix for k = 3 directly. As an 
easy consequence of the recursion formula, point ^ of Proposition 2.4, we get 
the following 

Corollary 6.1. The B-spline functions Ni = A^i,3, 1 < i < M of order 3 corre- 
sponding to the knot sequence A = (tj)|=i are given by 



(6.1) Ni{x) 



[X 



ti 



20),(10 



•2+2 



(20),(21), 



+ 



(X 



(31).(21), 



X] 



(31).(32), 







ifxe [ti,ti+i), 
ifxe [ti+i,ti+2), 

ifxe [ti+2,ti+3), 

otherwise. 
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13 



(6.2) 



u+2 /21)2 r 

N,{x)Ni+i{x)dx ^ 

ti+i 



(31). 



1 ^ (10). ^ (32), 



10 30(20)i 5(31)i 



and due to symmetry 



(6.3) 



iV,(x)iV,+i(x)dx=^' ' 



1 ^ (43). ^ (21), 



10 30(42)i 5(31) 
Proof. A straighforward calculation using Corollary |6.1 yields 

(21)2 r (^21) (10) 1 (21) 



Ni{x)Ni+i{x)dx 



20(20) 12(20) 4 5(31) 



An easy reformulation of this identity gives us (6.2). Equation (6.3) now follows 
by symmetry. 



□ 



Proposition 6.3. We have the following formulas for the Gram matrix A 

3- 

(32)? 



i^i,j)ij=i ~ ((-^i' -^i))i!j=i of B-splines of order 3. 



(6.4) ai,i+2 = ai+2,i 

(6.5) aj^j+i = cii+i,i 

(6.6) ai^i 



30(31),(42)i 

(31). , (21)f(10). ^ (32)f(43) 



10 30(20)i(31)i 30(31)^(42), 
(30), (30),(21)2 



5 15(20)i(31), 
Furthermore, = if \i — j\ > 2. 



if 1 < i < M - 2, 
if 1 < i < M - 1, 
ifl<i<M. 



Proof. Equation (6.4) is a simple consequence of Corollary 6.1 For equation (6.5), 
we observe that (A^,,A^,+i) is the sum of (6.2) and (6.3). Thus, (6.5) is a conse- 
quence of the identity 



(21)^ 
(31) 



1 _^ (32) 



10 5(31) 



(31) 



1 , (21) 
10 5(31) 



(21)2 + (32)2 (21)(32) (31) 



10(31) 



5(31) 



10 



where every bracket {mn) is taken with respect to the subindex i. We observe 
that from properties |3] and [5] of Proposition 2A we obtain 



i+2 



(6.7) 



N„ J2 Ni) = {Ni,l) 



(30)i 



l=i-2 



Since the only part in this equation that is unknown is (A^,, A^,), we use (6.7) and 
simple calculations to conclude the proof of the proposition. □ 

Observe that the Gram matrix A is symmetric, 2-banded and totally positive 
(see Remark 2.10), so in particular, we can apply Lemmas 4.2 and 4.4 to A. 
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6.2. Continuous piecewise quadratic splines. As a preliminary case, we ex- 
amine the setting 

(6.8) /i= (3, 2, 2,..., 2, 2, 3). 

where every knot in the mesh sequence is a double knot. Thus, the sequence 
{ti)^^i has the form 

(6.9) = ti = = ^3 < ^4 = ^5 < ^6 = ^7 < ■ ■ ■ • 



By Theorem |2.3[ this corresponds to the space of piecewise quadratic splines 
that are continuous at the mesh points Si,l < i < L. In fact, this is a special case 
of the choice fi = (3, 1, 1, . . . , 1, 1, 3) since the Gram matrix A of the B-splines 
depends continuously on the point sequence (tj)l=i (see for instance |Sch81j . The- 
orem 4.26). It follows from Proposition 6.3 and this continuity property that in 



the case of (6.8), formulas ( 6 . 4 ) - ( 6 . 6 ) for the Gram matrix A 
to 

(6.10) a,, 



■,jJi,j 



simplify 



(6.11) ai^i+i = {Ni,Ni+,) 

(6.12) ai^i+2 = {Ni,Ni+2) 




if i is odd 
if i is even 



if 1< z < M, 



10 

(32), 

30 



if 1 < i < M - 1, 
if 1< z < M - 2. 



There are two reasons why present the special case of (6.8). First, we obtain 
better estimates than the general ones of Section 6.3 Secondly, in contrast to 
Section 6.3, we do all calculations by hand. 



Remark 6.4. We note that (6.11) and (6.9) imply an^n+2 = if n is even. 
Furthermore, if n is even, (6.10) - (6.12) yield the following relations between the 
entries of A. 



(6.13) 



3a 



n— 1,71+1 



'^n,n+l) '^71—1,' 



''n,7i+l ; 



1; 3an,n — 4a 



n—l,n- 



Lemma 6.5. If A = {aij)fj^^ is defined as in (6.10) - (6.12) and 2 < n < M is 
even, we have 



(6.14) 



- -6" 

3 



for 1 < j < n. 



Proof. If n is even, the ra-th row of the matrix A has only three non-zero entries 
by Remark OA So the n-th row of An has only two non-zero entries and thus we 
have by the defining property of -B„ = {An)~^ 



n 



i=l 



/ , V^^i^i^n = an- l,n + &",nan,n for CVCry 1 < j < U. 

i=n—l 



The assertion of the lemma now follows from these equations and (6.13). 



□ 



Lemma 6.6. If A = {aij)fj^i is defined as in (6.10) - (6.12) and 2 < n < M is 
even, we have 



(6.15) 



6", 



LTL-l 

n— 1,71- 



-l"n-l,nJ 5 
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(6.16) 



_ 25 2 _ 14 



Proof. The first equation (6.15) is a consequence of (4.4) and tlie fact tliat tlie 



n-tfi row of An lias only two non-zero entries. For the second equation we observe 



that by (4.4) and the checkerboard property of Bn = (^n) ^ (see Remark 2.10) 



^n+l,n+l ~ {(^n+l,n+l b^,n'^n,n+l + '^\bn-l,n\^n,n+lO"n-l,n+l &J^-l,n-l'^n-l,n+l ) 



Using Lemma 6.5 for „ and b'^_^ we obtain 



^n+l,n+l ~ ( C^n+l.n+l - 



16 



^n,n+l ^(^n,n+lO'n-l,n+l + g Cln-l,n+l 

A 2 \ -1 



(^n—l,n 



3 a 



n— l,n 



By (6.13), this becomes 



tn+l 
"n+l,n+l 



''n+l,n+l 



and thus the lemma is proved. 



□ 



Proposition 6.7. If A = 

M, we have the estimate 



[a 



^ is defined as in (6.10) - (6.12) and 1 < n < 



(6.17) 

Proof. By induction. If n 



/ (10)n 



V 9 12 



(32) 

n 



(30), 



(10)i , (21)i , (32)i 



9 



12 



since (10)i = (21)i = 0. 

In the first place, we assume that n is even. In this case we get by equation 



(6.15) of Lemma 6.6 (6.10) -(6.12) and the inductive hypothesis respectively 

-.2 \-l 



un-l 2 N 

^n-l,n-l"n-l,ra/ 



^ / 2(21), 

V 15 
< /2(21), 



tn.-l 
"ri-l,n-l 



(10)„- 



< 



15 V 9 
/2(21)„ (21)„\-' 



V 15 



20 



(21). 
10 

^ (32)„_i 
5 
12 



-1 



Y(21)n^2 

V 10 



-1 



(21), 



/(10)„ (21)„ (32)„x-i 
V 9 12 ^ 5 y ' 



where we used that (10)„ = (32)„ = 0. Thus the claimed inequality (6.17) is 
proved for n even. 



Secondly, assume that n is odd. Here we get by equation (6.16) of Lemma 6.6 



(6.10) -(6.12) and the inductive hypothesis respectively 



6", 



_25 1 2 _14 N-i 

(^n,n g-j^ "n-l,n-l'^n-l,n 2Y^"'~^' 
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(30)„ 25^„_i ^(21)n-l^2 14(21)„_i 



< 



5 81 
(30)„ 25 12 



n— l,n— 1 



10 



27 10 



(21)„_a2 14(21)„_ 



5 81 (30)„_iV 10 ) 27 10 

(32>„x-i 



/ (30)n _ 4(10)„ \-i ^ /(lO)^ 

V 5 45 / V 9 

/(10)„ (21), (32)„x^i 



V 9 



where we used that (21)„ = 0. This proves inequahty (6.17) if n is odd. We thus 
proved the assertion of the proposition. 



□ 



Lemma 6.8. // A 

we have 



IS defined as m ( |(U0| ) - ( [(U2| ) and 2 < n < M, 



(6.18) 

if n is even and 
(6.19) 
if n is odd. 



Kn = -&n,n&L' i«n-i,„ for every 1 < j < n - 1 



6" 



"j,n-l"'n-l,n 



for every 1 < j < n — 1 



Proof. Equations (4.3) -(4.5) yield the formula 

in /tn— 1 



(6.20) 



Since for n even, a„_2,n = 0, equation (6.18) 



K,n{b],n-l^n-l,n + ^j,n-2'^n-2,n) • 

bllows. 



On the other hand, we get from Lemma 



6.5 



that — ^i,n-l/fln-2,n-l 



4&"'„ii/3 for 1 < j < n — 1. Combining this with (6.20), we obtain 



6" 

J," 



0'n-2,n-l 



Using (6.13) in this equation, we obtain (6.19). This completes the proof of the 
lemma. □ 



Proposition 6.9. If A = {aij)i^^-^^ is defined as in (6.10) -(6.12) and 2 < n < 
M , then the entries b'j^ of the last column of the matrix Bn = A~^ admit the 
geometric estimate 

(6.21) \Kn\<cS--^ for every 1 < j < n, 



r]{j,n) 

where q = (3/5)^/2 ^ ^ 



Proof. We can use Lemma 6.8, Proposition 6.7 and (6.11) to show that 



(6.22) 



ri{i,n) 



for n — 2 < j < n. 



where C,-,„ satisfies the inequality Cj^„ < 18 for all n — 2 < j < n. Indeed, for 



7 = n this is a direct consequence of Proposition |6.7[ If j = n — 1, we use Lemma 



6.8 



to write as an expression involving only 6^ ^, V^_\ and a„-i,n- Then we 



are able to invoke again Proposition 6.7 and equation (6.11) to show (6.22) for 
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j = n — 1. If j = n — 2 we apply Lemma |6.8|to write 6"^ as an expression involving 



only b'^ b^_2 „_i and fln-i.n- Here we use Proposition 6.7 for equation (6.11 ) 



for an-i,n and the already proved estimate (6.22) for 6J^_2„_i. 



If j < 72 — 3, we obtain from Lemma 6.8 and (6.11 ) 



Ife" I 



un-1 



hP- n n I 

^n-l,n-l^n,n^n-2,n-l^n-l,n\^j^n-2\ 



< 



5/ (10)„_i (21>„_ 
9V 9 

(31)n-2(31)n- l 

100 



12 



1 I (32). 



_A-i/(10)„ (21)„ (32)„x~i 
/V912 5/ 



< 



5 12-9(31)„_2(31)„_i , 
9 100 (30)„_i(30)„ ' ^■■"-2'' 



since at least one of the two terms (21)„_i, (21)„ vanishes. The elementary in- 

1 



equality (cf. Axiom (3) in Definition 2.5) 

(31)„_2(31),_i 



< 



(30)„_i(30)„r/(j,n-2) " 7]{j,n) 



induction and (6.22) now yield the assertion of the proposition. 



□ 



Given Proposition 6.9, we basically use the same proof as the one of Theorem 



5.3| to deduce the geometric decay off the diagonal for the inverse Gramian in the 
current case. This is formulated in the following 

Theorem 6.10. If A = {aij)f^j^^ is defined as m ( [6lo| - ( [6l2| ) andl<n<M, 
then the entries 6^- of the matrix Bn satisfy the geometric inequality 



(6.23) 



3C. 

vl/2 



for 1 < i, j < n, 



where, as in Proposition 6.9, q = (3/5) ' and C = 18 



Proof. Since Bn is symmetric and the case j = n was treated in Proposition [6^ it 
suffices to consider the case i<j<n — l. At first we only look at the parameter 
choices i < j and j odd. Equations (4.3) - ( 4^ yield the formula 

(6.24) bl^ = blj' + b-JlJbl^. 



Equation (6.24) and Lemma 6.8 imply with a„ = 1 if n is even and a„ = 5/9 if 
n is odd 



(6.25) 

By induction, 
(6.26) 



6". 



' "n^n,n^i,n-l"i,n-l"n-l,n- 



n-1 



i+i- 



We now use that ai < 1 for all / and employ Propositions 6.9 and 6.7 to deduce 



from (6.26) 



n-l 



l\(^l,l+l 
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Figure 1 : Definition of the Gram matrix A in M ATHEMATIC A. 



15 



(3-x[± + l]'/ ((x[i] +x[i + l]) (x[i + I] +x[i + 2]))) /; j == i; 



x[i + I] + x[i + 2] 
10 



a 



(x[i + I]^x[i]) / (30 (x[i] +x[i + l]) (x[i + l] + x[i + 2])) + 

(x[i + 2]^ x[i + 3] ) / (30 (x[i + I] + x[i + 2] ) (x[i + 2] + x[i + 3] ) ) // j == i + 1; 
j_] := x[i + 2]^ / (30 (x[i + 1] + x[i + 2] ) (x[i + 2] + x[i + 3] ) ) /; j == i + 2; 



Figure 2: Definition of the function ip in M ATHEMATIC A. 



<p[i_] : = 

1 1 1 1 

-x[i] + — x[i + I] + - x[i + 2] - (x[i + I] x[i + 2]) / (30 (x[i + I] + x[i + 2])) 
,5 12 5 

(1 + x[i + 2] / (5 (x[i + 1] +x[i + 2] ) ) - (20x[i] ) / (5 (x[i] + x[i + 1]) ) ) 

+ (5x[i- 1] x[i]) / (108 (x[i- 1] +x[i])) + 



(2 x[i - 1]' x[i]) / (73 (x[i-l] +x[i])^) 




r/(«,j) ^ 



n— 1 



12 .p, 



(31)? 



(30),+i r7(z,/) 77(j,/) 100 



C 



r/(2,j) V 25 0(30)m 



where q,C are as in Proposition 6.9 Since ri{i,j) < -^(z, /), (31)/ < ri{j,l) and 



(31); < (30)/+i for all j < / < n — 1, we get further 

(6.27) i*;y<^;^fi + ^;,« 



1 + ^ n2('-i) 



25 



U j = i or j is even, the proof follows the same lines, but — in view of Lemma 
6.8 — with a bit more notation, so we omit it. This is the place where the summand 
3 appears in (6.23) instead of 1 in the right hand side of (6.27). □ 



6.3. Continuously differentiable piecewise quadratic splines. In this sec- 
tion we come back to the general case for piecewise quadratic splines, where 
k = 3, fi = (3, 1, 1, . . . , 1, 1, 3) and the Gram matrix A of the splines Ni = A^j 3 



has the form (6.4) -(6.6). In the following, we use a computer algebra system, in 



our case Mathematica 8.0, to show that certain given polynomials have only 
nonnegative coefficients. This allows us to deduce that the given polynomial it- 
self is nonnegative for positive arguments. In the following results, we need the 
matrix A to be defined in Mathematica. This is done in Figure [T] Furthermore 



we need a Mathematica expression for to be defined in Proposition |6.11 
For this, we refer to Figure [2j 



The first thing to show is, as in Section 6.2, a suitable upper bound for b[ 



This is the content of the following Proposition 6.11 
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Proposition 6.11. Let A be as in Proposition 6.3\ and 1 < n < M. Then the 
bottommost rightmost element of Bn = A^^ satisfies the estimate 



(6.28) 



< ('(^Q) I (^^) I (3^) (21)(32) / (32) 20(10) x 
9 12 5 30(31) V 6(31) 9(20)7 



^ 5(0,-l)(10) ^ 2(0,-1)^(10) 



108(1,-1) ' 73(1,-1)2 
where every bracket (ij) in this formula is taken with respect to index n. 
Remark 6.12. Observe that in the hmiting case (10),^, (32)„ — )■ for n even, ipn 



is the estimate of Proposition 6.7 



Proof. We use Lemma 4^ and induction to prove the claimed estimate. Recall 
that Lemma 14.21 stated that 



"n.+ l,n+l — 



(6.29) 



Oin,n+lOin—l,n+l 



2(3^1, n'^n- 



71+1 



^n— l,n 



) 



Ojn—l,n 



"'?i-l,n+l'^n-l,n-ll-^ ^ ^n,n^n-l 



n— l^n— 1.71/ 



for 2 < n < M — 1. Thus to apply induction, we need to verify (6.28) for n = 1 
and n = 2. First, we suppose that n = 1; here we have = a^i = 5/(30)i by 

and some 



(6.6) and the fact that (20) i = 0. For n 



2, the formula b"^ 2 



«2,2 



holds in view of Corollary 4.1[ We thus obtain by (6.5), (6.6), (20)i 
elementary calculations 



(6.30) 



(21)2 ^ (32)2 (21)2(32)2 



12 



30(31)s 



(32)2 \ 
(6(31)^ 



The expression for b\ ^ and (6.30) are special cases of (6.28) since (20)i = 0. Thus 



the lemma is proved for n = 1 and n = 2. 

Before we proceed with the actual proof of (6.28), we show that the term 
ttnn+i — 2""."°"-i."+i ^ appearing in (|6.29l), is nonnegative. It is equivalent to show 

' ^n— l,n 1 1 

that an,n+iO'n-i,n — '^0'n,nCin-i,n+i is nouuegative. This is done using the Mathe- 
MATlCA-code of Figure |3| The method of proof is as follows. 

(1) Define the rational function p = an,n+ian-i,n — 2a„^„a„_i^„+i depending 
on the variables x[l] = (10)n_i,x[2] = (21)„_i, . . . ,x[5] = (54)„_i. 

(2) Determine the denominator d of the rational function p as 900 (x[l] + 
x[2])(x[2] + x[3])2(a;[3] + a;[4])2(x[4] + a;[5]) and observe that d is positive 
for positive arguments. 

(3) Calculate the coefficients of the polynomial q := d ■ p and verify that no 
coefficient of q is negative. 



Now we continue with the proof of (6.28). Since every entry a^j of the Gram 
matrix A is nonnegative and the matrices are checkerboard, the argu- 



ment of the last paragraph shows that a sufficient condition for (6.28) to be true 
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Figure 3: Proof of the inequality an,n+i — 2a„nan-i,n+i/an-i,n > 0. 

in[ii]:= p[z ] := Factor[a[l, 2] a[2, 3] - 2 a[2, 2] a[l, 3] I. 

3'ajble[x[n] -LLstyzWn^ , {n, I, Lengthl-LLstiz^^^W ; 
d[z ] := DenomxnatorypyzW ; 

g[z ] ■•= P[z] ; 

d[Seguence@@ Tajbie[x[n] , {n, 1, 5}]] 
Select [Flatten [Coefficientlist [ 

g[ Sequence @@ rafale [x[n] , {n, 1, 5}]], rai3le[x[n], {n, 1, 5} ] ] ] , B < fi] 

Out[U]= 500 (x[l] +x[2] ) (x[2] +x[3] )^ (x[3] +x[4] ) (x[4] +x[5] ) 

Out[15]= {} 



for all 1 < 72 < M is the following recursive inequality for 2<n<M — 1. 

/ 2ctfj ,2*^71— l,n+l \ 

(6.31) ■ ^ ^ 

(^n—l,n ^ 

The proof that this inequality is true for all choices of nonnegative distances 
(0, — !)„_!, (10)„„i, . . . , (54)„_i is equivalent to prove that the rational function 
p defined as 
(6.32) 

V^n-l'^n-l,n I (C^n+l.n+l ~ ^nO'n,n+l{0,n,n+l ) 

\ ^ Ojn—l,n 



and depending on the six variables (0, — l)„_i, (10)„_i, . . . , (54)„_i is nonnegative 
for nonnegative arguments. This is done in the MATHEMATlCA-code of Figure |4] 
using the following steps. 

(1) Define the rational function p to be the expression in (6.32) depending 
on the variables = (10)„_2, a;[2] = (21)„_2, . . . ,x[5] = (54)„_2, a;[6] = 
(65)„_2. 

(2) Determine the denominator d of the rational function p as an integer 
multiple of (a; [1] + a; [2] )^ (x [2] + X [3] f{x[^]+x[A]f{x[A] + x [5] ) ^ (a; [5] + a; [6] ) ^ 
and observe that d is positive for positive arguments. 

(3) Calculate the coefficients of the polynomial q := d ■ p and verify that no 
coefficient of q is negative. 

This proves the assertion of the proposition □ 

Remark 6.13. It is an easy consequence of the above lemma and the inequality 



^^^(^1 + ^U(32) 
(31) 1^ + 6(31)^-^^^^' 



that we also have the estimate 

< ^(10) ^ (21) ^ (32) 
(6.33) ^"-^^^ + T^ + ^J ='-'^- 

for all 1 < n < M. The Mathematica expression of t/'m is given by Figure [5] 
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Figure 4: Proof of the estimate for 

in[3i]:= p[z ] := Factor^ip[3]'^ ip[2]'^ a[2, 3] 

a[4, 4] -<p[3] a[3, 4] (a[3, 4] - (2a[3, 3] a[2, 4]) /a[2, 3]) - 

2a[3, 4] a[2, 4] ) 

a[2, 4]^<p[2] (l + ip[3] <p[2] a[2, 3]^) - <p[4]-' /. 

a[2, 3] ' j 

Table[x[n] List[z] ZnT\ , {n, 1, length [List [z] ]}]] ; 

d[z ] ;= Denominator [p[z'\ ] / 

g[z ] ■■= d[z] p[z] ; 

d [ Sequence @@ rafale [x[n] , {n, I, 5}]] 
Select [Flatten[CoefficientList [ 

g[Seguence@@ rafale[x[n] , {n, 1, 6}]], rafaie[x[n], {n, I, ff}]]],tt<0«] 

Out[34]= 72 441550 057 348 800 000 {x[l] +x[2])^ 

(x[2] +x[3])^ (x[3] +x[4])« (x[4] +x[5])^ (x[5] +x[5])^ 

Out(35]= { } 



Figure 5: Definition of the function ip. 



11 1 > 

— x[i] + — x[i + 1] + — x[i + 2] 
15 12 6 



Lemma 6.14. Let A be as in Proposition 6^ and 2 < n < M. Then we have the 
estimate 

6 (20)„ 



a 



n—l,n 



< 



5(30),' 



Proof. By Remark |6.13[ it suffices to show that 

6 (20)„ 



n~l.n 



< 



5 (30), 



We apply the same method of proof as in the above proposition and proceed with 
the M ATHEM ATIC A-code of Figure |6] using the following steps. 



1^1) Define the rational function p 



6 (20)^ 



5 (30)„ 



n— 1,71 



depending on the vari- 



ables x[l] = (10)„_i,x[2] = (21)„_i,x[3] = (32)„_i,x[4] = (43)„_i. 

(2) Determine the denominator d of the rational function p as 5(x[l] + x[2]) 
(x [2] + X [3] ) (x [3] + a; [4] ) (x [2] + X [3] + a; [4] ) (4x [2] + 3x [3] + 6x [4] ) and observe 
that d is positive for positive arguments. 

(3) Calculate the coefficients of the polynomial q := d ■ p and verify that no 
coefficient of q is negative. □ 



— ). This expression is important, since it is 
for 1 < j < n — 1. In the following lemma, 
product of two consecutive values of 9n- We will use this result 
in the proof of Proposition 6.16 



Let 9n := &"^„(a.n-i,n 
used in Lemma 
we estimate the 



4.4 



to estimate 



6" I 



to obtain explicit estimates for \W 
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Figure 6: Proof of the estimate for b^^an-i,n- 

^ 6 x[2] + x[3] 

in[2i]:= p[z ] := Factor i/f[2] a[l, 2] I . 

^ 5 x[2] + x[3] + x[4] 

rable[x[j] -» list[z] InJ, {n, 1, I,ength[List [z] ] } ] ] / 

d[z ] ;= DenomLnatorypYz^ ] / 

•JEz ] ■■= p[z] ; 

d[Seqiience@@ rable[x[n] , {n, 1, 4}]] 

SBlBCt\FlattenyCoBfficLent.'List.\_q\_x, y, z, w] , {x, y, z, w}]], tt < 4] 

Out[24]= 5 (x[I] +x[2] ) (x[2]+x[3]) (x[3]+x[4]) 

(x[2] + x[3] + x[4] ) (4 x[2] + 3x[3] + 5x[4] ) 

Out[251= {} 



Figure 7: Proof of the estimate for the product of two consecutive values of 



ln[261:= p[z ] 



Factor 



[- 

^ 100 



■<P[3] <p[4] 



x[3] + x[4] + x[5] x[4] + x[5] + x[6] 



a [2, 3] 



too 

a[I, 3] a[2, 2] 



x[3] + x[4] x[4] + x[5] 

a[2, 4] a[3, 3] 



a[3, 4] 



a[2, 3] 



/. 



a[l, 2] 

raJbIe[x[n] ^ List [z] InJ , {n, I, Length[L±st[z]]}]^; 

d[z ] := Denominator[p[z]] ; 

g[z ] ■■ = P[z] ; 

Select [Flatten [Coefficientiist [ 

d [ Sequence @@ rajble[x[n] , {n, I, ff} ] ] , rajb2e[x[n], {n, 1, ff}]]], B< Ofi] 
SeIect[FIatten[CoefficientI.ist[g[SeqTJence@@ ra£>Ie[x[n] , {n, 1, 6}]], 

raiile[x[n] , {n, 1, 6}]]] , n < &] 

Out[29]= { } 



Out[30]= { } 



Lemma 6.15. Let A be as in Proposition 6^ and 3 < n < M — 1. Then we have 
that 

87 (20)„ (20)„+i 



(6.34) 



100 (30)„(30)„+i- 
Proof. We show that the rational function defined as 

_8;^ _ (30)„ (30)„+i 
100 



(6.35) 



(20)„ (20)„+i 

depending on the variables x[l] = (10)„_2, a^[2] = (21)„_2, • • • , a;[6] = (65)n_2 is 
nonnegative for nonnegative arguments. This is done with the Mathematica- 
code in Figure [7] using the steps 

(1) Define the rational function p as in (6.35). 

(2) Determine the denominator d of the rational function p and verify that 
no coefficient of the polynomial d is negative. 

(3) Determine the coefficients of the polynomial q := d ■ p and verify that no 
coefficient of q is negative. □ 
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Proposition 6.16. Let A be as in Proposition 6.3, 1 < n < M and 1 < j < n. 

Then we have the estimate 



< c 



, , wi/i g = (87/100) C = 12g"^ 



6\2 



Proof. We see in the first place that Remark 6.13 



< 



12/(30)„ and thus the assertion of the theorem in t 



yields the estimate b" 
le case j = n. If j < n — 1, 



we get from Lemma 4.4 
(6.36) 



l=j+2 



Now assume that — (j + 2) + 1 = n — j — 1 is odd. By Lemmas 6.14 and 6.15 

12 /6\2 -n- (20), 



we conclude from (6.36) that 
(6.37) \b]J < 



(30), V 5 



1 



n-j~2 



n 

i=j+i 



mi 



We apply axiom (3) of Definition 2.5 on rj and induction to obtain the final 
estimate 

-5/ r]{j,n) 

If we assume that n — (j + 2) + 1 is even, the same proof yields the estimate 

_i6 g"--'' 



(6.38) 



\bln\<l2q-'Cl 



(6.39) 



IKnl < 12g- 



5r]{j,n)' 



Inequalities (6.38) and (6.39) together now prove the proposition. 



□ 



The passage to estimates of expressions for 1 < i,j < n is now an 



analogous calculation as in Theorem 5.3 and carried out in the following 



Theorem 6.17. Let A be as in Proposition 6.3 and 1 <n < M . Then the entries 
b^j of the matrix satisfy the estimate 



(6.40) 



for all 1 < i, j < n, 



where Ci = C(l + Hj^^-) and C, q are as in Proposition 



6.16 



Proof. Since Bn is symmetric and the case j = n was treated in Proposition 6.16 



it suffices to consider the cases i < j < n — 1. Equations (4.3) -(4.5) yield the 
formula 

" - t,j ' i,n j,n/ n,n 



By Lemma 4.4, inequality (4.10), we obtain 



\rri I ^ i 1;,"-^ h^-^ 



Applying Lemma 6.14 and Proposition 6.16 we estimate further 
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By induction, we get 



n-l 



(20) 



-,21-i-j 



1=3 



< C 



n— 1 
1=3 



(20) 



i+i q 



2(1-3) 



(30),+i '^'^vi^JHjJ) 



An easy consequence of formula (6.5) for a/^;+i is that < 2(31);/15 < 

2?7(i,/)/15. This and the obvious inequahties ri{i,j) < ri{i, I), (20) i^i < (30);+i 
for / as in the above sum give the final estimate 



b?A < C 



1,3 



■y3-^ 



12 
75 



C 



1=3 



1 + 



12 C 



751 



□ 



Observe that the Main Theorem 3.1 is a corollary to Theorem 6.17 Thus the 
proof of Theorem 3.1 is completed. 

7. Applications 

In the whole section, let A = {si)^^i be an arbitrary partition of the unit 
interval, > 2 an arbitrary integer and set ^ = (fc, 1, 1, . . . , 1, 1, k). Furthermore, 
o" = (A; — 1, . . . , — 1) and we denote by Pa the orthogonal projection operator 
from L^[0, 1] onto Sjt(A; a). Additionally we assume that the B-splines Ni = Ni^^ 
are constructed as in Section 12.21 and their Gram matrix and its inverse are 



denoted hy A = {aij)fj^^ and B = (&i,j)ij=i, respectively. Moreover, we define 
the Dirichlet kernel K corresponding to the space Sfc(A,(T) as 

M 

(7.1) i^(t,s) = ^6i,,A/,,(t)Ar,(s). 

For an arbitrary subset U of the unit interval [0, 1], we set A(f/) := A fl f/ 
understanding the partition A as a subset of the unit interval [0, 1]. Using A = 
(sj)fl]^, we define |A([/)| := max {sj+i — Sj:l<i<L — 1, {sj, Si+i} C f/}. 

Let / be an integrable function on the unit interval [0, 1]. We define the max- 
imal function M{f) of / as 



(7.2) 



M(/)(x) 



sup |/| 



\f{t)\dt for alia; e [0,1], 



and the maximal average function f* of / as 



(7.3) 



sup 



jj{t)&t for all X e [0,1], 



where in both cases, the supremum is taken over all subintervals / of [0, 1] con- 
taining X. Clearly, f*{x) < Mf{x) for all x G [0, 1]. 

A point a e [0, 1) is called a weak Lebesgue point of a locally integrable function 
/if 



(7.4) 



/(a) 



lim 



f{x) dx. 
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The Lebesgue differentiation theorem states in particular that a. e. point — with 
respect to Lebesgue measure — is a weak Lebesgue point. 

We remember the notation from Example 2.6 , where the function rji : {1, . . . , M} x 
{1, . . . , M} ^ M+ was defined by 

Viihj) = . max {k,0)i. 

iAj<l<t\/j 

We shall need the following partial integration lemma due to Taberski jTab62j . 

Lemma 7.1 ( |Tab62j ) . Let f G L^[a,b] and g be a function of finite total varia- 
tion on the interval [a + e, b] for each e > with 

(7.5) / vai{g; [s,b])ds < oo. 

J a 

Here and in the following, var(/; [a,b]) denotes the total variation of f on the 
interval [a,b]. 

Under these assumptions, the following limit exists and we have 
(7.6) 

\im f{t)g{t)dt < sup -/ /(t)dt / {\gib)\+v^Tig;[s,b]))ds. 

^^OJa+e 0<h<b~a H J a J a 



form 



The theorems in this section are stated for general k under the hypothesis that 



the inverse Gramian B = (6ij)f^=i satisfies a geometric decay inequality of the 



n\i-j\ 



for some constants < g < 1, Ci > 0. In Theorem 3A_ we proved such a bound 
with constants q, Ci independent of the underlying partition A for k = 3. This 
allows us to apply each result below for this choice of k, that is for piecewise 
quadratic splines. 

First, we present an ad hoc argument in the proof of Theorem 7.2 that ma- 
jorizes Pa/ pointwise by the maximal function M(/). This can be used to show 
that Pa/ converges back to / a. e. for every integrable function /, provided 
|A| ^ 0. 

Theorem 7.2. Suppose that 

(7.7) |&^,,|<Cli^ foreveryl<i,j<M, 

where < q < 1 and Ci > are some constants independent of i,j. Then we 
have that 

(7.8) \PAf{t)\<C2M{f){t) 

for all f e L^[0, 1] and t G [0, 1], where Ci = Ci ^ff is independent of f andt. 
Proof. Let t G [0, 1]. We first note that Pa is given by 

M 

(7.9) PA/(t) = ^6,,(/,iV,)iV,(t). 
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+ -t\)7]i{t,j) = + max {k,0)i 



(7.10) 



> ^ (^, 0); > t{iyj)+k — tiAj- 
l=iAj 

Now we choose the integer zq such that t G [tjo,tio+i) if t G [0,1) and we set 
io = M ii t = 1. Then we have by (7.9) and inequahty ( |7.7 ) 

jo iw 

I^A/(t)| = ' 



j=io — fc+l i=i 



A/ 



j=io— fc+l 



^i(«,J 



■|(/,iv,)l- 



Using (7.10) we see that 

to M 

(7.11) |PA/(t)|< £ Ar,(t)5^Ci 

j=io— fc+i i=i 



(|j -2| + l)gl*-J'l 



N,{s)\f{s)\ds. 



Since Ar,(s) < 1 pointwise and t G [tiAj,t{ivj)+k] if ^ and j are in the range of the 
sums in (7.11), this is less or equal 



M 



(7.12) 



J2 iV,(t)^Ci(|j-z| + l)gl'-^lM/(t). 

i=io — k+l j=l 



By Proposition 2.4, the B-splines Ni are a partition of unity, so it follows that 

M 



(7.13) 



|PA/(t)|< max VCidj 



z| + l)gl*~^'lM/(t). 



An elementary calculation gives + — ^(i^g)2 for all 1 < i < oo 



Consequently, we obtain the assertion of the theorem. 



□ 



Remark 7.3. If / G C[0, 1], Pa/ — ^ / uniformly as |A| — )■ 0. This is a conse- 
quence of the uniform boundedness of Pa as an operator from C[0, 1] to C[0, 1] for 
fixed spline order k (proved in |Sha01j ). Moreover, C[0, 1] is dense in L^[0, 1] and 
M is of weak type (1, 1). By Theorem 3.1, we may apply Theorem 7.2 to deduce 
|PA/(t)| < C2M{f){t) for A; = 3 and all / G L^[0, 1], where C2 is independent of 
A and /. From these facts, we obtain by a standard argument 

^lim^PA/(x) = /(a;), for a.e. a; G [0, 1], / G L^[0, 1] and k = ?>. 

In fact, this property is local in the sense that for every open subset U C [0, 1] 
we have 

lim PA/(a;) = fix), for a.e. x G ?7, / G LHO, 11 and = 3. 
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We will refine these results in the following in the sense that we bound Pa/ 
pointwise by the maximal average function /* and show convergence of -Pa/(o) 
to /(a) as |A| — > for all weak Lebesgue points a of /. 

The proofs of the theorems below are taken basically from |CK97] , where these 
results were shown for piecewise linear splines. 

The following theorem, concerning the total variation of the Dirichlet kernel, 
is essential. 

Theorem 7.4. Suppose that 

(7.14) \k,j\<Ci^-^^ for every l<z,j <M, 

where < q < 1 and Ci > are some constants independent ofi,j. Then, for all 
t G (0, 1), we have 



(7.15) max 



^ var(K(t,-);[s,l])ds,^ var(ir(t,-);[0,s])ds| <C3, 



where C3 = '^'^1i_qy?^ (md var(/; [a,b]) denotes the total variation of f on the 
interval [a,b]. 

Proof. We prove the inequality for the first expression in the maximum. The proof 
for the second expression is similar. We first fix t G [0, 1) and choose the integer 
io such that t G [tig,tig+i). Then we have by elementary properties of the total 
variation 

[ vai{K{t,-);[s,l])ds < [ var(ir(t, ■); [s, 1]) ds 

/ var(ir(t,-);[s,l])ds 

i=io 

rU+i 

<J2 var(ir(t,-);[t.,l])ds 

i=io 
M M 

<J2i^o),J2Y^i{K{t,-y,[t„tj+,]). 

i=io j=i 

We now observe that the total variation of B-spline functions Nj is < 2. This is 
a consequence of the following identity for derivatives of B-splines 

-(.W^(*-i)(^-|^), 



1=10 

(7.16) M 



since we have by this formula and Proposition 2.4 

var(Ar,,fc;suppAr,,fc) = / \Nl,,{t)\dt 

Jti 



{k 



u 

(A;,l), (A; -1,0), 



dt 
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k-1 



|^i+l,A;-l||Li(0,l) + 



k - 1 



j,fc-l||Li{0,l) 



We use this bound on the total variation of Ni^^ to estimate the total variation 
of the Dirichlet kernel appearing in (7.16) as 

H) j 

var(ir(t,-);[t„t,+i])<2 Yl E \bi,mmt) 

l=io—k+l m=j — fc+l 



The fact that the B-splines are a partition of unity and inequality (7.14) give us 



var 



(7.17) 



l=io—k+l m=j—k+l 

< 2k max \b^rn\ 

i()—k+l<l<io ' 
j—k+l<m<j 



< 2kCi max 



\m-l\ 



io-k+i<i<io rii{l, m) 

j—k+l<m<j 

<2kCiq-^^^q^-'° max r]i{l,m)-\ 

io—k+l<l<io 
j—k+l<m<j 



Observe that (10)j < m) for all /, m in the range of the above maximum and 
j > ^ > "^o- Thus we estimate the last line of (7.16) using (7.17) by 

M M 

(7.18) 

This concludes the proof of the theorem. 
Lemma 7.5. Suppose that 



□ 



(7.19) 



for every 1 < i,j < M, 



where < g < 1 and Ci > are some constants independent of Let a,s E 
[0,1) and define the integers io,jo such that a G [tjQ,tip+i) and s G [tj^, tj^+i). 
Then the Dirichlet kernel K satisfies the hound 



(7.20) 



\K{a,s)\<Ciq 



-fc+i g'^°'^°'(l + Ko-Jol) 
Is — a\ 



Proof. Since the B-splines are a partition of unity and due to estimate (7.19) we 
have 

M 



\K{a,s)\ = J2 KmNi{s)N„ 



(7.21) 



',m=l 



< max 

io — fc+l</<io 
jo-k+l<m<jo 



< Ci max 



\l-m.\ 



ia-k+l<l<io rii{l,m) 

jo-k+l<m<jo 
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Next we observe for /, m in the range of the maximum in (7.21 ): 

22) r]i{l,m) > _ max (10)„ > — ^ " 



«oAjo<"<«oVjo 



Ko - Jol + 1 



and 



\l — m\ > \io — jol — k + 1. 



This and (7.21) allow us to deduce 
(7.23) \K{a,s)\<Ciq 

This inequality proves the lemma. 
Theorem 7.6. Suppose that 



gKo-io|(l + 1 


'o - jol) 




S — 





□ 



(7.24) 



\Kj\<Ci 



for every 1 < i,j < M, 



where < g < 1 and Ci > are some constants independent of i,j. Then we 
have the inequality 

\PAfix)\<C,f*ix) 

forallf E L\0,1) and x E (0,1), where = 2{Cs + Ciq-''+^ maxi<,<M qK'^+'i')) 



is independent of f and x. Here, the constant C3 is the same as in Theorem \7.4 
Proof Let t E (0, 1) and / E L\0, 1). We first note that 



(7.25) 



PAfit) 



K{t,s)f{s)ds. 



Lemma 7.1 yields that 



K{t,s)f{s)ds 





1 


rt+h 


< sup 


-th 


Jt ^^'^ 


0<h<l~ 





{\K{t,l)\+YaTiK{t,-y,[s,l]))ds. 



We use now Theorem 7_A^ Lemma 7.5| and the definition of the maximal average 
function /* to deduce 



^ K{t,s)f{s)ds\ < {C,q 



-k+l 



m^x q\l + ^) + Cs)f*it). 



By a symmetric argument, we get the same inequality for the second integral 
K(t, s)f{s) ds in (7.25), and thus the theorem is proved. □ 

Theorem 7.7. Let f E L^{0, 1) and a be a weak Lebesgue point of f . Additionally, 
let U be an arbitrary neighbourhood of a. Suppose that 



(7.26) 



\h,\<Ci 



q' 



for every 1 < i,j < M, 



where < g < 1 and Ci > are some constants independent of i,j and A. Then 
we have 



(7.27) 
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Proof. We assume without loss of generality that U is an open interval containing 
a. We remark that K(t, s)ds = 1 for all t G [0, 1], thus we have to estimate 



(7.28) 



PAfia) - f{a) = C Kia, - /(a)) ds. 

Jo 



Let e > be an arbitrary positive real number. Since a is a weak Lebesgue point 
of /, there exists 5 > such that 



(7.29) M= sup 

0<\h\<5 



a+h 



fit) - f{a) dt <e and {s : \s - a\ < 5} (Z U. 



We now split the integral in (|7.28|) into the three parts 
(7.30) 



a+5 pa 

+ 1 + 

\s—a\>5 Ja J a—S 



We begin with the first integral Jig_^ 



(7.31) 



\s-a\>5' 

< ll/lli sup \K{a,s)\ + \f{a)\ [ \K{a,s)\ds. 

\s-a\>5 \s-a\>5 J\s-a\>S 



If \A{U)\ 0, the expression supj^^^j^^ |-ft'(a, s)| converges to zero by Lemma 
Now we estimate the integral J\g_a\ys l-^l*^' enough to consider the 



7.5 



integral jj^^ |i^'(a, s)| ds, since the other part follows by symmetric arguments. 
First we let the indices io, jo be such that a G [tio,tjo+i) and a + 6 E [tjo,tjQ^i). 
Then we conclude by ([7^ and ([7^ 



/ \K{a,s)\ds<J2 ma,s)\ds 



M 



< 



mi 



maxjo<„<i(10)„ 



-A;+l+|«-io| 



=J0 



The right hand side of the above inequality converges to zero, provided |A(f/)| — 
0, so we have proved that 



(7.32) 



\K{a,s)\ds^^ if |A(f/)| 0. 



|s-a|>(5 



In the second step, we estimate the integral ^'^^^ Kia,., s){f{3) — /(c^)) ds. The 
third integral f^_g Kia, s){f{s) — /(a))ds is treated in the same way. We apply 
Lemma |7.1| to this integral with the choice a = a, b = a + 6, f replaced by the 
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function /(■) — /(a) and g(t) = K{a,t). This yields 
(7.33) 



a+S 



K{a,s){f{s)-f{a))ds 



< sup 


1 




h 


I 


0<h<S 





a+h 



{fit) - f{a)) dt 



a+S 



\K{a,a + 6) \ + var(ir(a, ■); [t,a + 6])dt. 



By (|7g), supo<,<5 il fit) - f{a) dt\ < e. Theorem yields 

fa+5 



/ vax{K{a,-)-[t,a + 5])dt<C^ 

J a 



and Lemma 7.5 



gives 



Inserting these facts in ( |7.33[ ), we get 
(7.34) 



K{a,s){f{s)-f{a))ds 
for some constant C independent of A and /. 



< Ce 



Combining (7.32) with (7.34) in the splitting (7.30), we see that \PAf{a)—f{a)\ 
is arbitrarily small, provided |A(t/)| is chosen sufficiently small. This proves the 
theorem. □ 

Remark 7.8. As a last remark we note that Theorem 13.11 can also be used to 
prove unconditionality for orthogonal spline systems of order 3 with arbitrary 
knots in L^, 1 < p < oo. This will be investigated in a future paper. 
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